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Abstract 

For the first time accurate measurements of electron and positron fluxes in the 
energy range 0.2-^10 GeV have been performed with the Alpha Magnetic Spectrom- 
eter (AMS) at altitudes of 370-^390 km in the geographic latitude interval ±51.7°. 
We describe the observed under-cutoff lepton fluxes outside the region of the South 
Atlantic Anomaly (SAA). The separation in quasi-trapped, long lifetime (O(10 s)), 
and albedo, short lifetime (O(100 ms)), components is explained in terms of the 
drift shell populations observed by AMS. A significantly higher relative abundance 
of positrons with respect to electrons is seen in the quasi-trapped population. The 
flux maps as a function of the canonical adiabatic variables L, a Q are presented for 
the interval 0.95 < L < 3, 0°< a a < 90° for electrons (E<10 GeV) and positrons 
(E<3 GeV). The results are compared with existing data at lower energies. The 
properties of the observed under-cutoff particles are also investigated in terms of 
their residence times and geographical origin. 
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1 Introduction 



Evidence for high energy (up to few hundred MeV) electrons and positrons trapped below 
the Inner Van Allen Belts has been published during the last 20 years. The existing exper- 
imental data in the energy range of 0.04-=-200 MeV come from satellites covering a large 
range of adiabatic variables [|23|, 0, [JO, [II]]. Additional information, at relatively higher 



energies, is furnished by balloon-borne experiments however these data cover a 

more limited spatial range and have larger uncertanties due to the shorter exposure times 
and the presence of background from atmospheric showers. Although the magnetic trap- 
ping mechanism is well understood, a complete description of the phenomena, including 
the mechanisms responsable for the injection and depletion of the belts as well as those 
determining the energy spectra is lacking, particularly for energies above a few hundred 
MeV. At lower energies, models are available for leptons and protons [p2[|,P|] based on 
the data provided by satellite campaigns, which are continuosly updated for istance in 
the context of the Trapped Radiation ENvironment Development project [P^J]. 

At higher energies the existing data come from measurements carried out by the 
Moscow Engineering Physics Institute. These data, taken at altitudes ranging from 



300-^1000 km with different instruments placed on satellites and the Mir station [|24 
0,0], established the existence of O(100 MeV) trapped leptons both in the Inner Van 
Allen Belts (stably trapped) and in the region below (quasi-trapped), and determined their 
charge composition [||, ||]. At these altitudes, the shell structure is strongly distorted 
in the vicinity of the SAA, and consequently the observations are sensitive to different 
regions of trapped particles: the Inner Van Allen belts over the SAA and quasi-trapping 
belts outside of the SAA. An example of the shell structure relevant at these altitudes 
is shown in Fig.|I]: the shell evolves essentially above the atmosphere which it intercepts 
around the SAA. 

The Russian measurements concern mainly the region of the SAA; very little data is 
available at the corresponding altitudes outside the SAA. The measured ratio of e + to 
e~ is found to depend strongly on the observed population type. In the SAA, electrons 
dominate the positrons by a factor ~ 10, a ratio similar to that observed for the cosmic 
fluxes, while outside the SAA the two fluxes are similar and comparable to the e + flux 
inside the SAA [[§]]. However, the situation is not completely clear, since other groups 



report a lower e excess (~ 2) for the SAA [[13 



In the following, we use the high statistics data sample collected by the AMS exper- 
iment in 1998, for a detailed study of the under-cutoff lepton fluxes in the 0(1 GeV) 
energy range. The data are analyzed in terms of the canonical invariant coordinates char- 
acterizing the particle motion in the magnetic field: the L shell parameter, the equivalent 
magnetic equatorial radius of the shell, the equatorial pitch angle, ao, of the momentum 
p with the B field, and the mirror field B m at which the motion reflection occurs during 
bouncing [ fll7| , fQ]. 
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Figure 1: Example of geometrical surface of a drift shell in a quasi-trapping belt. Notice 
the typical not-closed structure, in the vicinity of the SAA (SPENVIS package |27|). 



2 AMS and the STS-91 flight 

The Alpha Magnetic Spectrometer (AMS), equipped with a double-sided silicon microstrip 
tracker, has an analyzing power of BD 2 =0.1A Tm 2 , where B is the magnetic field intensity 
and D the typical path length in the field. A plastic scintillator time-of-flight system 
measures the particle velocity and an Aerogel Threshold Cerenkov counter provides the 
discrimination between proton and e . In the present analysis, a fiducial cone with a 28° 
half-angle opening aperture was defined to select the leptons entering the detector, resul- 
ting in an average acceptance of ~160 cm 2 sr. Further details on the detector performance, 
lepton selection and background estimation can be found in || and references therein. 

The AMS was operated on the shuttle Discovery during a 10-day flight, beginning 
on June 2, 1998 (NASA mission STS-91). The detector, which was not magnetically 
stabilized, recorded data during 17, 6, 7, and 14 hours pointing respectively at 0°, 20°, 
45°, and 180° from the local zenith direction. The results presented here are obtained from 
the data of these periods. The orbital inclination was 51.7° in geographic coordinates, at 
a geodesic altitude of 370-=-390 Km. Trigger rates varied between 100 and 700 Hz. The 
data from the SAA is excluded in our analysis. 

The shuttle position and the AMS orientation in geographic coordinates were provided 
continuously during the flight by the telemetry data. The values of L, a and B m of the 
detected leptons were calculated using the UNILIB package [ |P7|j l with a realistic magnetic 
field model, including both the internal and the external contributions [[]i~5|, fIT| . 

The AMS Field of View (FoV) in the (L,ao) coordinate space is determined both by 
the orbit parameters (geographic locations and flying attitude) and the finite acceptance 
of the detector. 

A simulation was developped to determine the AMS FoV along the orbit and evaluate 
the effects due to the finite detector acceptance. The results are shown in Fig.0. 
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Figure 2: Comparison of the field of view of AMS with balloons and satellite measure- 
ments in (L,a G ). In the insert plot, the AMS coverage in f3 Q vs a Q is shown. 



The («o, L) coverage is similar for the different attitudes. The finite acceptance of 
the detector influences essentially the definition of the lower contour. The upper limit is 
imposed by the orbital altitude and is described by the relation sinao = ^Jo.311/L 3 B m , 
where B m = 0.225G is the minimum mirror field encountered along the AMS orbit. Since 
the particles which are mirroring above the AMS altitude cannot be observed, particles 
with large equatorial pitch angles can only be detected at low L values (L<1.2). At larger 
L, only particles with a smaller a can be observed. Because of the fixed flight attitudes, 

the azimuthal (3$ coverage in the local magnetic reference frame (z=B, x=(VB)±, y=zxx) 
was not complete, as shown in the insert plot of Fig.^. 



3 Data Analysis 

To reject the cosmic component of the measured lepton fluxes, the lepton trajectories 
in the Earth's magnetic field were traced using a 4 th order Runge Kutta method with 
adaptive step size. The equation of the motion was solved numerically and a particle was 
classified as trapped if its trajectory reached an altitude of 40 km [||], taken as the dense 
atmosphere limit where the total probability of interaction is 50%, before its detection in 
AMS. Although satisfactory in most cases, this approach is less stable when the particle 
rigidity falls in the penumbra region, close to the cutoff value. In this case, the trajectories 
become chaotic and small uncertainties in the reconstructed rigidity and in the B field can 
lead to a misclassification. The validity of the adiabatic approach requires the parameter 



e = p/R to be small [|14|], p8|], where p is the equatorial Larmor radius of a particle 
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1411 shows that the motion becomes 



and R the field radius of curvature at the equator, 
cahotic if e > 0.1. The AMS data are consistent with this limit even though the detected 
particle energies are relatively high. To avoid such effects, we have defined an effective 
cutoff, -R e //, as the maximum rigidity value at a given magnetic latitude 8 m for which 
no traced lepton was found to be of cosmic origin. The R e ff values as function of the 
magnetic latitude are shown as filled triangles in Fig. [| We rejected from our sample all 
particles with R>R e ff. 
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Figure 3: Effective cutoff as function of the magnetic latitude (filled triangles), defined 
as the lowest rigidity below which no primaries are found. The highest rigidities above 
which all the leptons are primaries (points) and 50% are secondaries (empty triangles) 
are also shown. 

The residence time, T/, of the under-cutoff particles is computed, i.e. the total time 
spent by each particle in its motion above the atmosphere, before and after detection. 
The geographical location where the trajectories intercept the atmosphere determine the 
lepton's production and impact points, defined as the position from which the particle 
leaves or enters the atmosphere. 

The residence time distribution as a function of energy is shown for positrons in Fig.|3]; 
the same behaviour is observed for electrons. All observed leptons have residence times 
below ~ 30 s, with 52% of the e~ and 38% of the e + having a Tj < 0.3 s independent of 
their energy. The corresponding impact/production points are spread, for both e + and 
e~, over two bands on either side of magnetic equator, as indicated by the the yellow 
bands in Fig.||. 

A scaling law, Tf ~ E -1 , is observed for the remaining leptons: they are disposed in 
two diagonal bands separated by a difference in T/ of ~ 2.2/ E s. The impact /production 
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Figure 4: Residence time vs energy for e + . The same structure is observed for e 
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Figure 5: Geographical positions of production and impact points in the atmosphere. 
Yellow bands show the distribution for short-lived e~ , red/blue bands show the produc- 
tion/impact distribution for long-lived e~. A similar but complementary structure is 
observed for e + . 
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Figure 6: Distribution of intersection points with atmosphere for the drift shells crossed 
by AMS. Yellow region corresponds to shells with B m > 0.48 L 041 G, blue region to 
B m < 0.48 L 41 G. 



points for e + are localized in the red/blue spots of Fig.^: the same regions describe 
respectively the production/impact regions of e~. 

The data have been previously published by the AMS collaboration [[Q]] using the 
terminology of short-lived and long-lived to classify the particles with Tf below and above 
0.2 s. However, no interpretation was advanced at that time to describe the observed 
distributions. |i~6| l has discussed the AMS results qualitatively. 



An exhaustive explanation must take into account the geometry of the shells relevant 
to the AMS measurements and the fact that these shells evolve partially under the atmo- 
sphere; therefore, no permanent trapping can occur. The residence times are determined 
by the periodicities of the drift (r^) and bouncing {jb <C Td) motions, the type of motion 
which dominates depends on the relative fraction of the shell mirror points lying above 
the atmosphere. 

The impact/production points correspond to the intersection of the shell surfaces with 
the atmosphere, as shown in Fig.[| where particles generated in interactions are injected 
into the shells. Long-lived and short-lived particles move along shells with different values 
of B m or, equivalently ao, which determine the mirror height on each field line. For 
high B m values, or low ao, the mirror height is very low and the shells penetrate into 
the atmosphere at nearly all longitudes. Therefore particles are absorbed shortly after 
injection in the shells. This is shown by the yellow bands in Fig.|] corresponding to 
shells with B m > 0.48 L 041 Gauss; they reproduce well the impact /production points for 
short-lived leptons. When B m is lower, or «o closer to 90°, the shells descend below the 
atmosphere only in the vicinity of the SAA, indicated by the blue regions in Fig.|3], which 
correspond to shells with B m < 0.48 L 0,41 Gauss, and reproduce the impact /production 
points of the long-lived leptons. These particles can drift nearly an entire revolution 
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Figure 7: Flux maps for two different energy bins: A), B) e + , e between 0.315<E<0.486 
GeV and C), D) e+ e" between 1.77<E<2.73 GeV. 

before absorption in the atmosphere. For the short-lived component, the bouncing motion 
is dominant; the residence time is given by Tf = kr b , where r b is the bouncing motion 
period, and k is an integer or half-integer between 1/2 < k < 5. In the dipolar 
field model, r& is given by r b = f(a )L, where / is a slowly varying function of a a ; the 
upper limit is r ~ 300 ms for the AMS data. The sub-structure seen in the short-lived 
component of Fig.f| is due to the discrete values of k and the different L shells crossed 
during the AMS orbits. For the long-lived component, the drift motion is dominant and 
Tf = k' Td, where Td is the drift motion period, Td = f'(a )/ EL, where /' is a slowly 
varying function of a Q , and k! is a number less than one, corresponding to the fraction 
of a complete drift shell spanned by a particle. The two bands seen for the long-lived 
component of Fig.^ correspond to fractions of ~0.65 and ~0.25 of a complete drift. 

4 AMS Results 

For the description of under-cutoff fluxes, the energy E, the L parameter and the equatorial 
pitch angle «o were used (this is preferred to B m since limited to 0° -j- 90°). A three- 
dimensional grid (E, L, «o) was defined to build flux maps. A linear binning in ao and 
logarithmic variable size for L and E bins were choosen to optimize the statistics in each 
bin. The interval limits and bin widths are listed in Table 1. 

The flux maps in (L, «o) at constant E give the distribution of particle populations at 
the altitude of AMS. Nine maps at constant E have been made. Two different maps for 
two different energy bins of e + and e~ are shown in Fig.^. 

The flux is limited by the cutoff rigidity R c : on a given shell only particles with R 
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Figure 8: The integrated flux maps for e + (A) and e (B) between 0.205<E<2.73 GeV. 
The line shows the curve below which no quasi-trapped leptons are found. 



< R c are allowed to populate the shell, hence lower and lower energy particles populate 
higher and higher shells. 

The e + , e~ flux maps and their ratio in the energy interval 0.2-=-2.7 GeV are shown 
in Fig. |8| and Fig. |9] respectively. The solid line in the two plots identifies the lower 
boundary in (L,a Q ) below which no leptons can be found with residence times larger than 
0.3 s and is defined by the relation sina c = 0.8L -17 . 

Above the curve, for increasing values of a a , the long-lived component of the fluxes 
becomes increasingly dominant. This is demonstrated in Fig. [l^ where the same distri- 
butions, integrated over a Q (C,D) and L (A,B), are shown. The contributions of leptons 
with Tf < 0.3 s and Tf > 0.3 s are represented with dashed and solid lines respectively. 
Above a Q > 60° the flux is due substantially to the long-lived component; the e + flux 
represents ~ 80% of the total leptonic flux, while at the same level or less than the e~ 
flux in the low a Q region. The long-lived component dominates only at very low L values 
where the positron excess is more pronounced. 

This is seen clearly in the energy spectra for particles with «o > 70°, shown in Fig. [IT], 
which is superimposed with the lower energy measurements from MARYA [||]. At large 
pitch angles, the e + flux is higher than e~ flux by a factor ~ 4.5, in contrast with MARYA 
data which indicate the same level of flux for both lepton charges. The critical pitch angle 
a c can explain the presence of the two well separated components in the residence time: 
the albedo (short-lived) and quasi-trapped (long-lived) ones. Particles inside a cone with 
a half-opening angle a c around B will enter the atmosphere every bounce and therefore 
will disappear rapidly, while those outside it will enter the atmosphere only near SAA. In 
this context, a c can be defined as the equatorial bouncing loss cone angle, i.e. the largest 
pitch angle for which particles enter the atmosphere every bounce. Taking into account 
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Figure 9: The integrated e + /e ratio between 0.205<E<2.73 GeV. The line shows the 
curve below which no quasi-trapped leptons are found. 







i i i i i i _ 




























10 
























1 






, i , ,: , i , - 







20 40 60 30 20 40 60 80 
«a (Peg) (Pgj) 



-1 1 


'C') : 




D) : 














1 






i 




1 


1.5 2 2.5 


31 
L 


1.5 2 2.5 


L 



Figure 10: The integrated flux as function of «o and as function of L for e~ (A, C) and 
e + (B, D) between 0.205<E<2.73 GeV. The full line shows the long-lived component, the 
dashed line shows the short-lived component, while the points show the total flux. 
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Figure 11: Energy spectrum comparison between AMS and Marya for e + and e for 
particles with «o > 70°. In the insert plot, the e + / e~ ratio comparison is shown. 



a c , the residence time can be written as Tf = krb9(—a + a c ) + k'Td0(a o — a c ), where 6 is 
the Heavyside step function. The two terms correspond to very different components of 
motion (r& <C r^), according to the bounce loss cone on that field line. Furthermore, all 
the observed particles are in a drift loss cone angle since all of them enter the atmosphere 
within one revolution after injection. This explains the absence of a peak at 90° in the 
pitch angle distributions of Figs. [Ll] A and B. 



5 Discussion 

The AMS data establish the existence of leptonic radiation belts, with particle energies 
in the range of several GeV, below the Inner Van Allen belts. The particles populating 
these belts are not stably trapped since the corresponding drift shells are not closed over 
above the atmosphere in the region of the SAA. 

At any given L, a critical value of the equatorial pitch angle, a c , the bouncing loss 
cone, can be defined to distinguish the long-lived, or quasi-trapped, and the short-lived, or 
albedo, components of the fluxes. The same value is found to separate the regions where 
the e + /e~ ratio is above or around unity: the charge composition shows a clear dominance 
of positively charged leptons in a definite region of the (L,a ) space above a c (L). 

The observed behaviour distinguishes these belts from the Inner Van Allen belts and 
limits the possible injection/loss mechanisms to those acting on a time scale much shorter 
than the typical particle residence times. Mechanisms related to Coulomb scattering, like 
pitch angle diffusion, are ruled out since they imply much longer time scales. Moreover, the 
observed charge ratio distribution provides an important constraint for potential models. 

The interaction of primary cosmic rays and inner radiation belt protons with at- 
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mospheric nuclei in the regions of shell intersection with the atmosphere are a natural 
mechanism for the production of secondary leptons through the tt — \x — e or 7r - 7 — e 
decay chains. This leads to a e + excess over e~ and seems suitable to explain the observed 
charge ratio for the quasi-trapped flux [p5|,fl0[|l. However, for the albedo flux the charge 
ratio is of the order of unity, as seen in Fig. |^, and other mechanisms might be present. 

Recent Monte Carlo studies based on this mechanism have been able to fully reproduce 
the under-cutoff proton spectrum reported by AMS [jS|], while a less good agreement for 
the under-cutoff lepton spectrum was obtained. In fllBH the influence of geomagnetic 
effects, mainly related to the East- West asymmetry for cosmic protons, is taken into 
account to qualitatively explain the observed charge ratio. However, more refined studies 
are needed to definitely exclude contributions from other mechanisms, i.e. acceleration 
processes acting on the leptons resulting from the decays of /3-active secondary nuclei and 
neutrons of albedo and solar origin [p6|]. 

In conclusion, the AMS under-cutoff lepton spectrum can be described naturally in 
terms of the canonical adiabatic variables associated with the Earth's magnetic field taking 
into account the role played by the atmosphere. There are clear indications that ir decays 
can account for the quasi-trapped component of the flux, while the situation is less clear 
for the albedo component where other processes may contribute. 
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